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ABSTRACT 

We describe and report progress on the joint U.S. and German campaign 
to map the X-ray emission from the entire Cygnus Loop supernova remnant 
with the ROSAT High Resolution Imager. The Cygnus Loop is the prototype 
for a supernova remnant that is dominated by interactions with the interstellar 



medium and supplies fundamental physical information on this basic mechanism 
for shaping the interstellar medium. The global view that these high- resolution 

CO 
O , 

fx- . interstellar medium and the pivotal nature of cloud-blast wave interactions 



in determining the X-ray morphology of the supernova remnant. While 
investigating the details of the evolution of the blast wave, we also describe the 
interstellar medium in the vicinity of the Cygnus Loop, which the progenitor 
star has processed. Although we do not expect the X-ray observations to be 
complete until September 1997, the incomplete data combined with deep Ha 
images provide definitive evidence that the Cygnus Loop was formed by an 



explosion within a preexisting cavity. 
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1. Introduction 

Supernova explosions are one of the most basic processes that define the physical 
and chemical state of the interstellar medium (ISM). The interaction between supernova 
remnants (SNRs) and the ISM is complex and symbiotic. X-ray observations are important 
in the study of SNRs because thermal emission from hot material (kT ~ 0.1-1 keV) is 
the primary tracer of gas that the supernova remnant blast wave heats and accelerates. 
High resolution X-ray observations are important for studying SNR gas dynamics for 
three reasons. First, the length scale, £, for cooling and recombination behind a radiative 
shock is small: £ ~ 2.5 x 10 17 t^ 7 /n cm, where n is the preshock number density, and 
v s7 is the shock velocity in units of 100 km s _1 (McKee 1987). Second, the preshock 
medium is not uniform, having clouds and substructure ranging over scales from tens of 
pc and downward. Third, gas-dynamical instabilities are important, and multi-dimensional 
hydrodynamic calculations predict that fully developed turbulence occurs when the blast 
wave collides with an interstellar cloud (Klein, McKee, & Colella 1994; Klein, McKee, & 
Woods 1995). Consequently, structure is present down to the smallest scales that can be 
resolved computationally. The improved ROSAT-HRI sensitivity and spatial resolution 
compared to Einstein make it invaluable for exploring cloud-blast wave interactions. The 
ROSAT-HRI combines high 0.1-2.4 keV sensitivity (approximately three times better than 
the Einstein HRI; Zombeck et al. 1990), high spatial resolution (FWHM « 6" on axis), 
and high contrast imaging (the half-energy radius is 3" on axis, with scattering less than a 
few percent; Aschenbach 1988) to provide a powerful tool for investigating cloud-blast wave 
interactions. 

The region of the Milky Way towards Cygnus is viewed along the Carina- Cygnus 
arm and exhibits a complex array of bright and dark clouds on the Palomar Sky Survey. 
Most remarkable are the arcs and loops of faint Ha nebulosity, the smallest and brightest 
of which is the Cygnus Loop (£ = 74.0, b = —8.6). The Cygnus Loop is the prototype 
middle-aged SNR and is a particularly suitable candidate for a case study being nearby 
(770 pc; Minkowski 1958) and bright with low foreground extinction [E(B — V) = 0.1], 
thus ensuring unrivaled signal-to-noise and spatial resolution. From an observational 
perspective the Cygnus Loop is unique in terms of X-ray brightness, shock front resolution 
(1" e 1.1 x 10 16 cm), and interaction with a range of interstellar conditions. The Cygnus 
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Loop is also at sufficiently high galactic latitude that confusion with background emission 
from the plane of the galaxy is minimized. Unlike other well-known younger SNRs such as 
Cas A, Tycho, and Kepler, the structure of the Cygnus Loop is dominated by its interaction 
with its environment. The large diameter of the Cygnus Loop (3° = 40 pc) assures that it 
is interacting with different components of the ISM, from molecular gas on the northwest 
edge (Scoville et al. 1977), to diffuse atomic gas along the northeast and eastern limbs 
(DeNoyer 1975; Hester, Raymond, & Blair 1994), and low density, hot, ionized gas to the 
south (Ku et al. 1984). Previous studies of the Cygnus Loop have focused on individual 
regions selected for their brightness or morphological peculiarity (e.g. Hester & Cox 1986, 
Hester et al. 1994, Fesen, Kwitter & Downes 1992, Graham et al. 1995, and Levenson et 
al. 1996). The goal of this high- resolution investigation (Graham & Aschenbach 1996) is to 
make an unbiased survey and thereby establish globally the state of the medium into which 
the blast wave expands. A coherent picture emerges from our preliminary results (§2). We 
present a unified interpretation of the environment and development of the Cygnus Loop 
(§3), constrain the global models that are relevant (§4), and summarize the conclusions in 
§5. 



2. Observational Strategy and Data Reduction 

Our primary goal is to map the entire Cygnus Loop remnant at the highest spatial 
resolution that can be attained with the ROSAT-HRI. Achieving this ambition is not simple, 
because the spatial resolution is limited by the intrinsic instrumental properties and photon 
counting statistics. The on- axis resolution is about 6", but the off-axis imaging performance 
is dominated by telescope aberrations, which increase with increasing field angle. At 20', 
the edge of the useful field of view of the HRI, the resolution has degraded to about 30". 
The large range of surface brightness allows us the luxury of mapping the brightest regions 
at full HRI resolution, but forces us to seek a compromise between effective resolution and 
observing time in the dimmest regions. To give a concrete example, the surface brightness 
observed with the ROSAT HRI yields a count rate of 3 x 10 -5 — 1 x 10~ 3 s _1 in a 5" x 5" 
pixel. For the brightest regions along the shell (within a factor of two of the highest surface 
brightness) a 15 ks exposure will net 15 counts and 0.4 background counts. It would take 
longer than 450 ks to obtain a similar number of counts in the faint interior. We therefore 
adopted the following observing strategy. We use the full 6" resolution to observe the 
brightest regions, which restricts the effective field of view to 8' radius. Each of these 
24 pointings requires 15-20 ks to obtain 15 counts per pixel. In a comparable exposure 
time, the moderate surface brightness interior regions yield about 10 counts per 12" x 12" 
pixel. At this resolution, the effective HRI field of view has a 12' radius. We have observed 
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these nine fields for 20-25 ks. Observing the low surface brightness regions, we obtain 10 
counts per cell only when using a 15" or larger pixels, for which the full HRI field of view 
is suitable. Each of these 10 fields requires a 35 ks exposure. Finally, one very low surface 
brightness region along the western edge of the "breakout" requires a 50 ks exposure to 
obtain reliable counting statistics. The locations of the HRI pointings are shown on the 
Einstein IPC map (Fig. 1). 

A few fields were observed during the initial operation of ROSAT, and a concerted 
effort to make a complete map began during A05, in October 1994. Because of the large 
total observing time this program requires, it is a collaborative project between U.S. and 
German investigators. So far 771 ks out of a total 1066 ks of data have been obtained and 
we expect completion in A07, by September 1997. The complete mosaic will consist of 
47, fields of which 37 have been observed at least in part. We present processed data for 
26 fields, of which 15 are complete. The target names, field centers, and exposure times 
are listed in Table 1, and the observations that are presented here appear first. Column 4 
contains the live time, which is the useful portion of the observations, and column 5 lists 
the remainder of the requested time. 

The data were processed and combined into a single image using software designed to 
analyze extended sources and the diffuse background (Snowden et al. 1994). Over such 
large spatial scales, each observation is a distinct projection of the sky onto the plane of 
the detector, and these must be re-mapped to a common projection. For each observation, 
using only the accepted time intervals listed in the raw data files, three images are created. 
The first map contains the total counts, the second map is a model of particle background, 
and the third is an exposure map, which includes the variations in the detector quantum 
efficiency These images have 5" pixels, matching the instrument's angular resolution. 
These maps from individual fields are summed on a pixel-by-pixel basis into three mosaics 
of the entire Cygnus Loop field. For the initial analysis we present here, we use a uniform 
pixel scale and field of view from each observation. Only the central 17' of each observation 
in order to avoid the significant degradation of the point spread function and decreased 
efficiency at large off-axis angles, except for two fields where a radius of 20' was used to 
avoid a gap in the resulting image. The count rate map (Fig. 2) has 15" pixels and is the 
difference of the total counts and the total background, divided by the net exposure (Fig. 
3). Background subtraction is not a problem for the bright and moderately bright fields, 
but the background does become comparable with the source surface brightness for the 
lower surface brightness regions. The extent of the observations is demonstrated in Figure 
4, where the individual image boundaries are drawn on the count rate map. 

To display the high resolution information contained in the mosaic, we also show three 
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smaller fields that have been binned with 6" pixels (Figs. 5a, 6a, and 7a). For comparison, 
we present aligned optical images of the same regions taken through a narrow Ha filter 
(Figs. 5b, 6b, and 7b). These were obtained using the Prime Focus Corrector on the 0.8 m 
telescope at McDonald Observatory and are described in Levenson et al. (1997). After basic 
reduction consisting of bias subtraction and flat-fielding, individual images were re-mapped 
and averaged using the same technique that was applied to the X-ray observations. Figures 
5c, 6c, and 7c are color versions of these combined data sets to show the relative spatial 
distributions directly. 

3. Results 

The X-ray morphology of a SNR is a primary diagnostic of the medium into which the 
blast wave propagates. At first sight the X-ray morphology depicted in the ROSAT-HRI 
mosaic presents a bewilderingly complex view. This shows prima facie that the Cygnus 
Loop is dominated by its interaction with an inhomogeneous environment. In contrast, we 
note also that the Cygnus Loop appears nearly circular except for the breakout region to 
the south, which suggests that the SNR expands in a uniform medium. Any model for the 
Cygnus Loop must resolve this apparent paradox. 

There is considerable large- and small-scale azimuthal variation around the perimeter 
of the Cygnus Loop and significant structure projected across its face. X-ray emission is 
bright near the rim, but the brightest peaks all lie interior to the circle that circumscribes 
the X-ray edge. The two most prominent X-ray regions, NGC 6960 on the western edge and 
NGC 6992 to the northeast, lie at the largest projected distance inside this circle. These 
regions contain complex networks of X-ray filaments and high surface brightness clumps. 
The emission here is strongly limb brightened. Figures 8 and 9 show radial cuts of the 
emission from these regions. The limb-brightening (i.e. the ratio of brightness at center to 
edge) is about a factor of 6 in the northeast (Fig. 8) and 30 in the west (Fig. 9). The 
adiabatic expansion of a supernova blast wave in a uniform medium results in some degree 
of limb brightened X-ray emission but is inadequate to account for these profiles. We have 
calculated the X-ray surface brightness profile (Fig. 10) using the Sedov density profile 
(Sedov 1993), taking X-ray volume emissivity to be proportional to density squared, and 
assuming that it is independent of temperature, which is reasonable for the energy coverage 
of our observations (Hamilton, Sarazin, & Chevalier 1983). The maximum limb-brightening 
is 3.3 and occurs near the edge of the SNR. The clear message from this comparison is 
that the observed X-ray profiles in these bright regions of the Cygnus Loop have peaks 
in the wrong place and too great in contrast to be the result of an adiabatic blast wave 
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propagating in a uniform medium. Thermal evaporation does not explain the X-ray profiles 
either, as the evaporative solutions of White & Long (1991) demonstrate (Fig. 10; see 
also Cowie, McKee, & Ostriker 1981). The effect of thermal conduction is to provide more 
mass farther behind the blast wave, so these solutions always exhibit less extreme limb 
brightening than the adiabatic case. The Cygnus Loop is neither a Sedov remnant nor a 
Sedov remnant modified by thermal evaporation. 

3.1. X-ray and Optical Correlations 

Optical line emission alone reveals some details of shock structure. For example, faint 
Ha filaments trace the progress of the blast wave in neutral gas, [O III] emission reveals 
where postshock gas first starts to cool and recombine, and [S II] and bright Ha together 
show where complete cooling and recombination zones have developed. Utilizing X-ray 
and optical information together, however, provides a more complete understanding of the 
structure and development of supernova remnant shocks. 

The observed correlation between X-ray and optical line emission falls into three broad 
categories that are distinguished by surface brightness and morphology: (1) regions of 
high surface brightness X-ray emission distributed in long coherent arcs that are aligned 
tangentially with, and lie interior to, the bright optical filaments; (2) localized indentations 
in the shell, typically less than 20' in size, where bright X-ray emission occurs at large 
distances behind the projected edge of the blast wave; and (3) regions of low surface 
brightness, limb-brightened X-rays bounded on the outside by faint, narrow Balmer-line 
filaments. 

Examples of type (1) morphology are NGC 6992 in the northeast (a = 20 h 56 m , 5 = +32°, 
J2000) and NGC 6960 on the western edge (20 h 46 m , +31°). Type (2) regions are found 
on the eastern edge (20 h 57 m , +31°), the southeast cloud (20 h 56 m , +31°), and the bright 
funnel-shaped section (20 h 50 m , +32°) that is located about 30' northeast of the optical 
emission region known as the "carrot" (Fesen, Blair, & Kirshner 1982). These regions tend 
to be more turbulent and include the most compact and highest surface brightness knots 
within the Cygnus Loop. Categories (1) and (2) occur where the shock has encountered 
a greater than average column of gas, which has caused substantial deceleration of the 
blast wave (Hester & Cox 1986; Hester et al. 1994; Graham et al. 1995; Levenson et 
al. 1996). The difference between regions (1) and (2) may be due to viewing geometry, 
the shape of the cloud which is being encountered, and the time since enhanced density 
was encountered by the shock. The western edge provides a vivid example of a blast 
wave-cloud interaction because this region is viewed almost exactly edge-on, minimizing the 
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complications of projection effects (Levenson et al. 1996). The best examples of type (3) 
filaments occur at the extreme north (20 h 54 m , +32°), to the south of the southeast knot, 
and to the north of NGC 6960. These filaments show the location of the current projected 
edge of the blast wave because Balmer-line filaments occur where previously undisturbed 
atomic gas is shocked (Chevalier & Raymond 1978; Chevalier, Kirshner, & Raymond 1980; 
Hester, Danielson, & Raymond 1986). The X-rays do not completely overlap the optical 
emission but are offset to the interior. We measure offsets of 5 — 20" between the exterior 
edge of X-ray emission and the edge of the Balmer filaments; offsets of 10" (equivalent to 
10 17 cm) are typical. The separation between the X-ray and optical emission is likely due 
to the time required to heat the gas behind the shock front before X-rays are emitted. 
For a shock of velocity v s = 400 kms" 1 , an offset of 20" implies a heating time of about 
200 yr. This is consistent with Coulomb collisions alone equilibrating the postshock ion 
and electron temperatures (Spitzer 1962). If the blast wave has decelerated in dense gas 
where the Balmer filaments are detected, the same separation from the X-ray edge implies 
that a longer time has elapsed, so Coulomb equilibration alone is still a possible heating 
mechanism. Whether Coulomb collisions or some other mechanism heat the post-shock 
electrons, long equilibration timescales are expected, as Laming et al. (1996) demonstrate 
in analysis of a non-radiative shock observed in SN 1006. 



4. A Global View 

4.1. Cavity Explosion 

Various X-ray studies have successfully discussed individual regions of the Cygnus Loop 
in terms of blast wave-cloud collisions (e.g. Hester & Cox 1986, Hester et al. 1994, Graham 
et al. 1995, Levenson et al. 1996), and it is well established that large-scale collisions 
fundamentally determine the morphology of portions of this supernova remnant. However, 
the examination of particular regions does not lead to an understanding of the global 
environment in which the SNR evolves. Specifically, it is important to determine whether 
the clouds that have been struck by the blast wave are isolated and part of a random 
distribution or part of a coherent structure, such as an H II region cavity or wind blown 
bubble formed by the supernova progenitor. Previous global models of the Cygnus Loop 
(e.g. Falle & Garlick 1982; Tenorio-Tagle, Rozyczka, & Yorke 1985) ignore the evidence for 
and effects of encounters with clouds and instead model a SNR breaking out of a molecular 
cloud. On the contrary, these X-ray observations emphasize that over large scales the blast 
wave is running into dense material, not expanding into a contrasting lower density region. 

With the data presented here, it is clear that the Cygnus Loop is the result of a cavity 
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explosion. The observations of the interior are consistent with blast wave propagation 
through a rather smooth, low-density medium, and at the periphery, the edge of the cavity 
has been encountered. The enhancements in the ROSAT X-ray data show that the blast 
wave is currently interacting with denser material — denser than the medium through which 
it previously propagated — over 80% of the perimeter that we have observed. Some of these 
enhancements are narrow, tangential filaments, while others are the large structures that 
had been identified in earlier observations. Such extensive clumps of material cannot be 
characteristic of the pre-supernova ISM within the current remnant and thus distinguish 
the cavity boundary. We observe the cavity walls in these data! Only two extended sites, 
near a = 20 h 48 m , 5 = +32°10', and a = 20 h 57 m , 5 = +30°30', show no evidence for blast 
wave interaction with the cavity wall. 

The cavity walls are nearly complete, consisting of a neutral shell as well as large clouds. 
When the blast wave reaches the dense atomic shell, it decelerates. The Balmer-dominated 
filaments directly trace the motion of the blast wave in the atomic gas of the shell. The 
X-ray signature of this category (3) morphology (cf. §3.1) is enhancement at the SNR edge, 
where gas is piled up behind the slowing shock. Figures 8, 9, and 11 provide examples 
of this characteristic X-ray profile in three different locations, where the enhancement 
occurs at radii of 0.98i? snr , 0.98i? snr , and 0.96-R snr , respectively. These enhancements are 
not due to the global expansion of the blast wave in a uniform medium or evaporation 
in a uniform distribution of clouds. They are too slight, too limited in radial extent, and 
misplaced relative to the Sedov and evaporative model predictions. The Balmer filaments 
and associated X-rays are detected over half the perimeter of the SNR, indicating that blast 
wave deceleration in atomic material is common. Some regions exhibit the type (3) X-ray 
signature but the corresponding optical filaments are too faint to be detected. Although 
the blast wave edge is not directly observed in these cases, the prevalence of category (3) 
and the observed X-ray morphology imply that the blast wave is decelerating in the dense 
cavity wall in these locations, as well. 

A radial profile across the northwestern edge (20 h 46 m , +31°40') emphasizes the type (3) 
X-ray limb-brightening and avoids confusion with emission due to interactions with large 
clumps (Fig. 11). There is no extended bright optical emission due to radiative shocks near 
the edge, and this is not a particularly bright region of the X-ray map. Two components 
are obvious, corresponding to two overlapping projections of the blast wave that are not 
necessarily physically related. We neglect the interior component and model the exterior 
component very simply as a shock that strikes a density discontinuity and generates a 
reflected shock. The Rankine-Hugoniot jump conditions determine the relationship among 
the flow variables in different regions. The free parameters of the model are the density 
contrast and its initial location, from which the density profile due the projected and 



- 9- 



reflected shocks as a function of radius is calculated. In the Cygnus Loop, Balmer-dominated 
filaments are typically observed propagating through regions of n ~ 1 cm~ 3 (e.g. Hester 
et al. 1994, Raymond et al. 1983), while in the former H II region cavity n ~ 0.1 cm -3 , 
so the density contrast is taken to be 10. (This is likely to be an overestimate, however, 
because the brightest Balmer filaments were analyzed, and these are the ones that have 
evolved most toward the radiative stage.) The radius of the SNR blast wave is fixed at the 
location of the Balmer emission. The interior edge of the X-ray enhancement marks the 
location of the reflected shock. We measure this to be at R = 0.94i? snr , so the blast wave 
initially encountered the wall at R = 0.96R snr and is now at R = 1.0R snr . Similar to all the 
model calculations presented here, the X-ray surface brightness profile is computed for the 
spherically symmetric case (Fig. 11), assuming that X-ray volume emissivity is proportional 
to density squared and is independent of temperature. This simple model reasonably 
matches the shape and location of the X-ray enhancement, although this is not a fit to the 
data. The model does not allow for heating time, so it predicts emission toward the edge 
of the SNR blast wave that we do not detect. Despite these approximations, however, the 
simple model indicates the scales of the parameters that cause the X-ray enhancement at 
the edge where the blast wave has run into the atomic cavity wall. 

Several observed properties of the Cygnus Loop argue for a large-scale cavity explosion. 
The first is the approximate spherical symmetry of the emission. This symmetry of the 
Cygnus Loop is surprising in view of the rich X-ray morphology, which argues for an 
inhomogeneous medium. If the inhomogeneity were in the form of a random distribution of 
clouds, then the clouds would have to be very small in order to account for the fact that 
the observed gas produces few, if any, significant indentations into the remnant (McKee & 
Cowie 1975). There are no entirely concave sections on scales greater than 20' (4 pc) and 
no evidence for thermal evaporation of large clouds on the interior, for example. However, 
X-ray studies of individual regions show that the clouds are quite extended along the 
periphery of the remnant, with sizes of about 0?5 (7 pc). The clouds are preferentially 
aligned parallel to the blast wave. In this case, a spherical remnant can be produced only 
if neither the shape nor the distribution of the clouds is random. To produce a spherical 
remnant, the clouds in the vicinity of the Cygnus Loop must be located at the edge of the 
current shell, and they must be flattened along the shell. There is also dynamical evidence 
that these are the first large clouds that have been encountered by the blast wave: the 
blast wave-cloud interactions to the northeast and the west are recent and approximately 
the same age (Hester et al. 1994, Levenson et al. 1996). This conclusion is consistent with 
that of Charles, Kahn, & McKee (1985), who determined that data are inconsistent with an 
explosion in a uniform distribution of clouds based on Einstein-IPC observations adjacent 
to NGC 6960. They suggested that the ISM within the current filamentary ring lacked 



-10- 



clouds before the supernova explosion. 

4.2. Progenitor Modification of the ISM 

The modification of the SNR environment by the pre-supernova star has been discussed 
by McKee, Van Buren, & Lazareff (1984) and Shull et al. (1985). If the progenitor star 
of the Cygnus Loop SNR were an early B star, then the Lyman continuum radiation from 
this main-sequence star created an H II region. In the warm neutral interstellar medium 
with average density distribution n m 0.2 exp(— z/220 pc) cm -3 (Kulkarni & Heiles 1988), 
a star at height z = 120 pc emitting Lyman continuum photons at a rate of 10 46 s _1 
creates an H II region with a density of about 0.1 cm" 3 and a radius of about 30 pc. As 
denser clumps of gas within the H II region are photoevaporated, adding more material to 
the ionized interior, the radius of the H II region decreases. Small clouds are ultimately 
destroyed by photoevaporation and larger clouds are driven away by the rocket effect 
(Bertoldi & McKee 1990). The details of the photoevaporation depend on the velocity of 
the progenitor star and on the properties of the clouds: temperature, density, and magnetic 
field. Using Bertoldi's (1989) results, one can show that atomic clouds (T ~ 10 2 K) in 
pressure balance with the H II region (P/k ~ 2000 K cm" 3 for n# ii ~ 0.1 cm" 3 ) will be 
either destroyed or rocketed out of the H II region. Molecular clouds are more difficult to 
destroy by photoevaporation (Shull et al. 1985), but at the distance of the Cygnus Loop 
below the plane, very large molecular clouds are unlikely. As a result of photoevaporation 
and photodissociation, a dense atomic shell that contains most of the mass surrounds the 
H II region, which is bounded by a recombination front. A vigorous wind may aid formation 
of the shell, although this is unlikely in the case of an early B star we consider here. Once 
the progenitor leaves the main-sequence and becomes a red giant, the H II region begins 
to recombine, but the density is too low for much recombination to occur before the star 
explodes. 

When the progenitor of the Cygnus Loop became a supernova, the blast wave initially 
expanded adiabatically through the diffuse cavity gas. Eventually it hit the walls of the 
cavity, which are comprised of the surrounding atomic shell and the remains of clouds 
that protrude into the H II region. NGC 6960 is a clear example where the blast wave 
has encountered dense molecular gas (Levenson et al. 1996). The eastern edge and the 
southeast cloud may be additional examples of elephant-trunk structures, such as those 
observed in M 16 (e.g. Hester et al. 1996), protruding into the shell (Graham et al. 1995). 
The existence of large clouds near the cavity edge is expected, for the development of an 
H II region tends to homogenize its interior, while its action on the boundary is limited by 
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intervening clouds along the same line of sight and decreased ionizing flux at larger radii. 
The clouds currently at the edge were either rocketed out of the interior or were originally 
located at the periphery of the H II region. 

We have considered a progenitor of spectral type later than BO because a stationary 
star of earlier spectral type star would homogenize an H II region to a radius of about 50 
pc, while there are clearly inhomogeneities within the present 20 pc extent of the SNR 
(McKee et al. 1984; Charles et al. 1985). If a massive star is moving, however, the clouds 
situated perpendicular to its path are not exposed to the Lyman continuum emission for 
as long. The rocket effect and photoevaporation are less effective on these clouds, and 
the homogenizing radius across the direction of motion is reduced. Along the direction 
of motion, the distribution of the processed interstellar medium is anisotropic: toward 
the original location of the progenitor, gas has been exposed to ionizing radiation, so it 
is more homogeneous and less dense, while the processing is limited in the direction that 
the star moves, so the ISM is clumpier and less dense there. The progenitor of the Cygnus 
Loop would have been later than BO, provided that it was not moving too rapidly; (t>* ^ 
a few kms _1 so that surrounding clouds are exposed to the ionizing radiation for a time 
R/v* ^ 6 x 10 6 yr). Alternatively, an earlier-type progenitor is possible provided that it was 
moving predominantly along our line of sight. In this case, the projected edge of the cavity 
would still be circular, yet have a smaller radius than one carved by a stationary star. 

5. Conclusions 

The density inhomogeneities around the Cygnus Loop determine not only the 
appearance of the SNR at optical and X-ray wavelengths, but also they fundamentally alter 
the blast wave development. The consequences are significant beyond this single example. 
Global models of the interstellar medium that assume the clouds around supernova 
remnants are insignificant or rely on statistical measures such as the number-diameter 
relation must be revised. 

The high-resolution, global view of the Cygnus Loop these data provide resolves the 
apparent paradox of the Cygnus Loop: a near-circular shape and evidence for interaction 
with large-scale inhomogeneities in the interstellar medium. The X-ray data make clear 
that the blast wave has recently encountered density enhancements, which we identify 
as the walls of a cavity. The H II region of the progenitor star worked to homogenize a 
spherical cavity with ionizing radiation, and the blast wave initially propagated through 
this cavity, maintaining its spherical symmetry. Marking the edge of the cavity is a dense 
shell of atomic gas and large clumps of the processed surrounding medium or formerly 
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interior clouds that were rocketed out. Portions of the blast wave have been decelerated in 
the smooth atomic shell, which Balmer-dominated filaments and slightly limb-brightened 
X-ray profiles trace. Other sections of the blast wave have swept up even larger column 
densities, resulting in very bright X-rays and optical emission. The observed features of the 
Cygnus Loop are expected rather than surprising in the context of stellar evolution that 
preceded the supernova. 

We thank the USRSDC for providing the ESAS software used to reduce the X-ray 
data. The optical images were obtained using the Prime Focus Corrector on the 0.8-m 
telescope at McDonald Observatory in collaboration with L. Keller and M. J. Richter. The 
research of CFM is supported in part by NSF grant AST95-30480. 
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Fig. 1. — ROSAT-HRI pointings superposed on the Einstein IPC image of the Cygnus Loop 
(Seward 1990). Solid circles indicate fields that are presented in this paper; dashed circles 
show unprocessed or proposed observations. 

Fig. 2. — ROSAT-HRI count rate map of the Cygnus Loop. This image is a mosaic of 41 
pointings in 26 positions. Each observation has been background-subtracted and exposure- 
corrected before being combined into this map, which has 15" pixels. It is scaled linearly 
from to 0.15 counts s -1 arcmin~ 2 . 

Fig. 3. — Exposure map of combined ROSAT-HRI observations, including detector 
variations. 

Fig. 4. — The boundaries of individual pointings are shown on the ROSAT-HRI count rate 
mosaic of the Cygnus Loop. The sharp, nearly circular edge of the supernova remnant is 
real, not an artifact of the limited extent of the observations. 

Fig. 5. — The northeast rim of the Cygnus Loop at 6" resolution. (a^HRI data, which exhibits 
three distinct types of X-ray emission: long filaments of bright emission, spots of extremely 
bright emission accompanied by indentations in the X-ray shell, and fainter limb-brightened 
filaments. (b)The northeast rim of the Cygnus Loop in Ha, scaled linearly. The optical 
image contains extended radiative filaments, small spots of bright emission, and long, faint 
filaments due to Balmer-dominated shocks, corresponding to the classes of X-ray emission. 
^False-color composite of HRI (green) and Ha (red) observations, illustrating the exact 
relation between the X-ray and optical emission. The optical image has been median-filtered 
to remove stars. 

Fig. 6. — The southeast region of the Cygnus Loop at 6" resolution. (a)RRI, (b) Ha, and 
(c) false-color composite, as in Figure 5. 

Fig. 7. — The northwest edge of the Cygnus Loop at 6" resolution. fajHRI, (b) Ha, and (c) 
false-color composite, as in Figure 5. 

Fig. 8. — HRI surface brightness as a function of radius across the bright northeastern rim 
of the Cygnus Loop. This radial profile is an azimuthal average over 2'. Very bright emission 
appears projected to the interior of the SNR and extends over more than 12% of its radius. 
While there is no significant X-ray emission at the edge of the SNR, fainter X-ray limb- 
brightening appears to the interior, at R = 0.98R snr . (Balmer-dominated optical filaments 
define the edge of the SNR that is indicated in Figures 8, 9, and 11.) This X-ray emission 
is clearly distinct from that due to expansion of a blast wave in a uniform medium (cf. Fig. 
10). 
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Fig. 9. — HRI surface brightness as a function of radius across the bright western rim of the 
Cygnus Loop, azimuthally averaged over 2'. Extremely bright emission appears projected to 
the interior of the SNR, at R = 0.92R snr . Significant, though less extreme, limb-brightening 
is also observed at R = 0.98R snr . 

Fig. 10. — Model X-ray surface brightness profiles for Sedov SNR (solid line) and evaporative 
cloud model of White & Long (1991; dotted line) with their parameters C = 10 and r = 10. 
Both profiles have been normalized to a uniform level at the center of the supernova remnant. 

Fig. 11. — HRI surface brightness as a function of radius across the northwestern rim of the 
Cygnus Loop, azimuthally averaged over 2'. Two separate components of the limb-brightened 
shell are evident, near R = 0.96R snr and R = 0.85R snr . Both components are the result of 
interaction of the blast wave with the cavity wall, although the latter appears projected to 
the interior of the SNR. The outer peak is modeled as the result of a density contrast and 
reflected shock (dotted line). In this model, the blast encounters a density contrast of 10 
at R = 0.96R snr and the reflected shock has progressed to R — 0.94R snr when the forward 
shock has reached R = 1.0R. nr 
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